Theory of nonlinear optical properties of phenyl-substituted polyacetylenes 
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In this paper we present a theoretical study of the third-order nonlinear optical properties of 
poly(diphenyl)polyacetylene (PDPA) pertaining to the third-harmonic-generation (THG) process. 
We study the aforesaid process in PDPA's using both the independent electron Hiickel model, as 
well as correlated-electron Pariser-Parr-Pople (P-P-P) model. The P-P-P model based calculations 
were performed using various configuration interaction (CI) methods such as the the multi-reference- 
singles-doubles CI (MRSDCI), and the quadruples-CI (QCI) methods, and the both longitudinal 
and the transverse components of third-order susceptibilities were computed. The Hiickel model 
calculations were performed on oligo-PDPA's containing up to fifty repeat units, while correlated 
calculations were performed for oligomers containing up to ten unit cells. At all levels of theory, 
the material exhibits highly anisotropic nonlinear optical response, in keeping with its structural 
anisotropy. We argue that the aforesaid anisotropy can be divided over two natural energy scales: 
(a) the low-energy response is predominantly longitudinal and is qualitatively similar to that of 
polyenes, while (b) the high-energy response is mainly transverse, and is qualitatively similar to 
that of trans-stilbene. 

PACS numbers: 78.30.Jw,78.20.Bh,42.65-k 
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I. INTRODUCTION 

Conjugated polymers are believed to be one of the 
most important materials for the future nonlinear opto- 
electronic devices^. The nonlinear optical response of 
these materials originates mainly from the ir electrons 
which (a) are delocalized along the backbone of the 
polymer (b) while being completely localized as far 
as the motion in a direction transverse to the back- 
bone is concerned. It is this one-electron picture of 
localized-delocalized ir electrons which provides an in- 
tuitive understanding of the large response of these ma- 
terials to external fields. However, the role of electron- 
correlation effects on the linear and nonlinear nonlinear 
optical response of conjugated polymers has attracted 
considerable amount of theoretical attention in recent 
years, and a number of conjugated polymers such as 
irans-polyacetylene, poly-(para)phenylene (PPP), poly- 
(para)phenylenevinylene (PPV) etc. have been stud- 
ied using the Coulomb correlated models in recent 
works ^iiSiSiLL2i From several of these studies the pic- 
ture that emerges is that the nonlinear optical properties 
of various conjugated polymers for the most part are de- 
termined by a small number of excited states. For conju- 
gated polymers whose structures are invariant under the 
inversion operation (centro-symmetric) these states are: 
(a) 1B U , the lowest-energy one-photon state, (b) mA g , a 
two-photon state higher in energy than IB U and with a 
strong dipole coupling to it, and (c) nB u , a one-photon 
state higher in energy than mA g (but still within the 
conduction-band threshold), and with strong dipole cou- 
pling to it&2Jl&ijL The question is whether this so-called 
"essential-states picture" is peculiar to a few conjugated 
polymers, or is it universally shared by all of them? 

Recently, optical properties of a new class of 
conjugated polymers called phenyl-disubstituted poly- 
acetylenes (PDPA's) — which are obtained by substitut- 



ing the side H atoms of irens-polyacetylene by phenyl 
derivatives — have been the subject of several experi- 
mental studiesi 10 ) 11 ) 12 ) 13 ) 14 ) 15 ) 16 ) 17 The interesting aspect 
of the linear optical properties of these materials is 
that they exhibit strong photoluminescence (PL), de- 
spite their structural similarities to irens-polyacetylene, 
which is well-known to be nonphotoluminescenti 1 ^ In 
a series of theoretical studies, we explained the strong 
PL exhibited by PDPA's in terms of reduced electron- 
correlation effects caused by the delocalization of exci- 
ton in the transverse direction because of the presence 
of phenyl rings in that direction i 1 ^ 1 !?! 2 ^ We also showed 
that delocalization of exciton also leads to reduction in 
the optical gaps of these materials as compared to trans- 
polyacetyleneii&AiLSfi Moreover, we predicted that the de- 
localization of excitons will also leave its signatures in 
form of a significant presence of transverse polarization 
in the photons emitted during the PL in PDPA's, a pre- 
diction, which since then, has been verified in oriented 
thin-film based experiments. 17 

In the present work, we have undertaken a system- 
atic theoretical study of nonlinear optical properties of 
PDPA's corresponding to the third-harmonic generation 
(THG) process. There are several motivations behind 
the present study, the first of which is studying those ex- 
cited states of PDPA's which will not be visible in the 
linear optics. In the THG spectrum of centrosymmet- 
ric materials such as PDPA's, A g -type states appear as 
two-photon resonances, while the B u -type states appear 
as three-photon resonances, thus making the simultane- 
ous exploration of both the symmetry manifolds possi- 
ble. The next motivation is to explore whether one can 
understand the nonlinear-optical properties of complex 
polymers such as PDPA's, in simple terms as enunci- 
ated in the essential-states mechanism mentioned above. 
In particular, for irans-polyacetylene, THG spectrum is 
predicted to have two prominent peaks corresponding to 
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the essential states mA g and nB u mentioned above, in 
addition to the 1B U peaki^ Therefore, questions arise: 
(a) How does the nonlinear optical response of PDPA's 
compare to that of irans-polyacetylene? (b) What is the 
nature of the essential states contributing to the non- 
linear optical properties of PDPA's? (c) Will transverse 
polarization also play a role in the nonlinear optical prop- 
erties as it did in case of their linear optical properties? 
The last, but not the least of our motivations behind the 
present study is to investigate whether one can use con- 
ventional computational methods such as configuration- 
interaction (CI) approach also to study large unit cell 
materials such as PDPA's (fourteen electrons per cell). 
This is important in light of the fact that with so many 
novel and increasingly complex materials being discov- 
ered by "nano-technologists" on a routine basis, there 
will be novel challenges for anyone interested in obtain- 
ing their theoretical understanding. Therefore, it is im- 
portant to apply time-tested theoretical approaches such 
as the configurations-interaction (CI) method to study 
these materials to explore whether they can withstand 
the challenges posed by systems of increasing complex- 
ity. We address all these issues by performing calcula- 
tions on oligo-PDPA's of the third-order nonlinear optical 
susceptibilities corresponding to the THG process using 
both the independent-particle Hiickel model, as well as 
Coulomb-correlated Pariser-Parr-Pople (P-P-P) model. 
The P-P-P model-based calculations were performed us- 
ing the CI methodology, and approaches such as multi- 
reference-singles-doubles CI (MRSDCI), and quadruples- 
CI (QCI) were used to compute the THG spectra. For 
the Hiickel model calculations, we considered oligomers 
of PDPA's containing up to fifty unit cells, while for 
P-P-P model calculations up to ten unit cell oligomers 
were considered. We computed susceptibility component 
where all photons involved were polarized in the chain 
direction (longitudinal component), as well as the com- 
ponent where the photons were polarized perpendicular 
to the chain (transverse component). We see that the 
THG spectrum is clearly divided over two distinct en- 
ergy scales: (a) the low-energy response is primarily con- 
centrated in the longitudinal component and is qualita- 
tively similar to that of irans-polyacetylene, while (b) the 
high-energy response is concentrated mainly in the trans- 
verse component and is similar to that of irans-stilbene. 
Upon examining the orbitals involved in these responses 
we conclude that it is the low-energy chain-like molecular 
orbitals (MOs) close to the Fermi level which are active in 
the longitudinal component, while the transitions among 
the chain-like and high-energy phenyl-based orbitals give 
rise to the transverse component of the response. 

The remainder of this paper is organized as follows. In 
section [H] we briefly describe the theoretical methodol- 
ogy used to perform the calculations in the present work. 
Next in section IIIII we present and discuss the calcu- 
lated nonlinear optical susceptibilities of oligo-PDPA's. 
Finally, in section IWl we summarize our conclusions. 



II. METHODOLOGY 

The unit cell of PDPA oligomers considered in this 
work is presented in Fig. ^ Ground state geometry 
of PDPA's, to the best of our knowledge, is still un- 
known. However, from a chemical point of view, it is 
intuitively clear that the steric hindrance would cause 
a rotation of the side phenyl rings so that they would 
no longer be coplanar with the polyene backbone of the 
polymer. The extent of this rotation is also unknown, 
however, it is clear that the angle of rotation has to 
be less than 90 degrees because that would effectively 
make the corresponding hopping element zero, implying 
a virtual disconnection of the side phenyl rings from the 
backbone. In our previous worksi&ISiSfi, we argued that 
the steric hindrance effects can be taken into account by 
assuming that the phenyl rings of the unit cell are ro- 
tated with respect to the y-axis by 30 degrees in such 
a manner that the oligomers still have inversion symme- 
try. It is obvious that along the direction of conjugation 
(longitudinal direction), PDPA is structurally similar to 
irans-polyacetylene, with alternating single and double 
bonds. But, perpendicular to the conjugation direction 
(transverse direction), the material has features in com- 
mon with irens-stilbene, with its two phenyl rings, and 
vinylene-like linkage connecting them. In the following, 
we will adopt the notation PDPA-n to denote a PDPA 
oligomer containing n unit cells of the type depicted in 
Fig. □ 

The point group symmetry associated with trans- 
polyacetylene (and polyenes) is C^u so that the one- 
photon states belong to the irreducible representation 
(irrep) B u , while the ground state and the two-photon ex- 
cited states belong to the irrep A g . Because of the phenyl 
group rotation mentioned above, the point group symme- 
try of PDPA's is Ci so that its ground state and the two- 
photon excited states belong to the irrep A g , while the 
one-photon excited states belong to the irrep A u . How- 
ever, to facilitate direct comparison with polyenes, we 
will refer to the one-/two-photon states of PDPA's also 
as B u /A g -type states. 

The independent-electron calculations on the 
oligomers PDPA-n were performed using the Hiickel 
model Hamiltonian which, adopting a notation identical 
to our previous worksi&i&SSi reads, 

H = Hc + H p + Hcp, (la) 

where He and Hp are the one-electron Hamiltonians 
for the carbon atoms located on the irans-polyacetylene 
backbone (chain), and the phenyl groups, respectively, 
Hcp is the one-electron hopping between the chain and 
the phenyl units. The individual terms can now be writ- 
ten as, 

H C = - {h-{-l) M ^t)B Kk ,, M ,M+u (lb) 

(k,k'),M 
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Hp — —to Bw,m,m, (lc) 

{H,v),M 

and 

Hcp = —t± ^ Bk^r.Ai.M,- (Id) 

{k,/j,),M 

In the equation above, k, k! are carbon atoms on the 
polyene backbone, /x, ^ are carbon atoms located on the 
phenyl groups, M is a unit consisting of a phenyl group 
and a polyene carbon, (...) implies nearest neighbors, 

and l!,, r .\i.M- = Y!,A c i,M,a c 3,M',<r + h.c). Matrix ele- 
ments t , and t± depict one-electron hops. In He, At is 
the bond alternation parameter arising due to electron- 
phonon coupling. In Hcp, the sum over fi is restricted 
to atoms of the phenyl groups that are directly bonded 
to backbone carbon atoms. There is a strong possibil- 
ity that due to the closeness of the phenyl rings in the 
adjacent unit cells, there will be nonzero hopping be- 
tween them, giving rise to a term Hpp in the Hamilto- 
nian above. However, in our earlier studyjA& we explored 
the influence of this coupling on the linear optics of these 
materials, and found it to have insignificant influence. 
Therefore, in the present study, we believe that we are 
justified in ignoring the phenyl-phenyl coupling. 

As far as the values of the hopping matrix elements are 
concerned, we took to = —2.4 eV, while it is imperative 
to take a smaller value for t±, because of the twist in 
the corresponding bond owing to the steric hindrance 
mentioned above. We concluded that for a phenyl group 
rotation of 30 degrees, the maximum possible value of t± 
can be -1.4 eV4& Bond alternation parameter At = 0.45 
eV was chosen so that the backbone corresponds to trans- 
polyacetylene with the optical gap of 1.8 eV in the long 
chain limit. 

The THG process in oligo-PDPA's were studied 
by computing the third-order nonlinear susceptibilities 
X^{—3ui;ui,ui,uj). For short we will refer to this suscep- 

(3) 

tibility, as Xthg- Fi fS t the Hiickel Hamiltonian for the 
corresponding oligo-PDPA was diagonalized to compute 
the one-electron eigenvalues and eigenfunctions. These 
quantities were subsequently used in the formulas derived 
by Yu and Supi to compute Xthg- 

The correlated calculations on the oligomers PDPA-n 
were performed using the P-P-P model Hamiltonian 

H = Hc + H P + H C p + H ee , (2) 

where He, Hp, Hcp have been explained above. H ee 
depicts the electron-electron repulsion and can be written 

as 

H ee = U n^Uji 

i 

+^E^>*- 1 )(^- 1 ) ( 3 ) 



where i and j represent all the atoms of the oligomer. 
The Coulomb interactions are parameterized according 
to the Ohno relationship^!, 

V i j=U/K i j(l + 0.m7F% J ) 1 '* , (4) 

where, k»,j-,m,jv depicts the dielectric constant of the 
system which can simulate the effects of screening, U is 
the on-site repulsion term, and Rij is the distance in 
A between the ith carbon and the the jth carbon. The 
Ohno parameterization initially was carried out for small 
molecules, and, therefore, it is possible that the Coulomb 
parameters for the polymeric samples could be some- 
what smaller due to interchain screening effects* 2 - 2 Since 
the results obtained will clearly depend on the choice of 
the Coulomb parameters, we tried two sets: (a) "stan- 
dard parameters" with U = 11.13 eV and Kij = 1.0, 
and (b) "screened parameters" with U = 8.0 eV and 
Ki,i,M,M — 1 -0, and Ki,j t M,N — 2.0, otherwise. Using 
the screened parameters, Chandross and Mazumdar 22 
obtained better agreement with experiments on excita- 
tion energies of PPV oligomers, as compared to the stan- 
dard parameters. Recently, we performed a large-scale 
correlated study of singlet and triplet excited states in 
oligo-PPV's and observed a similar trendi 2 ^ 

As far as the hopping matrix elements for PDPA's and 
polyenes are concerned, they were assigned the values 
used in the Hiickel-model calculations, except for the 
bond-alternation parameter for which the smaller value 
At = 0.168 eV, consistent with the P-P-P model, was 
used. In trans-stilbene, for phenyl rings hopping value of 
-2.4 eV was used, while in the vinylene linkage the values 
used were -2.2 eV for the single bond, and -2.6 eV for the 
double bond. In all the calculations, C-C bond length 
of 1.4 A was used for the phenyl rings. In polyenes 
and PDPA's, along the backbone the single bonds and 
the double bonds were taken to be 1.45 A and 1.35 
A, respectively. The bond connecting the backbone to 
the substituent phenyl rings was taken to be 1.40 A. In 
trans-stilbene, in the vinylene linkage, the single (dou- 
ble) bond lengths were taken to be 1.54 A (1.33 A). The 
bond lengths used in this paper are the same as the ones 
used in our earlier works ^i^Q^ 

The starting point of the correlated calculations for 
various oligomers were the restricted Hartree-Fock (HF) 
calculations, using the P-P-P Hamiltonian. The many- 
body effects beyond HF were computed using differ- 
ent levels of the configuration interaction (CI) method, 
namely, quadruples-CI (QCI), and the multi- reference 
singles-doubles CI (MRSDCI). Since the number of elec- 
trons in oligo-PDPA's is quite large despite the P-P-P 
approximation owing to the large unit cell, it is not pos- 
sible to include all the orbitals in the many-body cal- 
culations. Therefore, one has to reduce the number of 
degrees of freedom by removing some orbitals from the 
many-body calculations. In order to achieve that, for 
each oligomer we first decided as to which occupied and 
the virtual orbitals will be active in the many-body cal- 
culations based upon: (a) their single-particle HF en- 
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ergies with respect to the location of the Fermi level, 
and (b) Mulliken populations of various orbitals with re- 
spect to the chain/phenylene-based atoms. Because of 
the particle-hole symmetry in the problem, the numbers 
of active occupied and virtual orbitals were taken to be 
identical to each other, with the occupied and virtual or- 
bitals being particle-hole symmetric. The remaining oc- 
cupied orbitals were removed from the many-body calcu- 
lations by the act of "freezing", i.e., by summing up their 
interactions with the active electrons, and adding this 
effective potential to the one-electron part of the total 
Hamiltonian. The inactive virtual orbitals were simply 
deleted from the list of orbitals. When we present the CI 
results on various oligo-PDPA's, we will also identify the 
list of active orbitals. During the CI calculations, full use 
of the spin and the point group (C, for PDPA's) symme- 
tries was made. From the CI calculations, we obtain the 
eigenfunctions and eigenvalues corresponding to the cor- 
related ground and excited states of various oligomers. 
Using the many-body wave functions, we compute the 
matrix elements of the dipole operator amongst various 
states. Finally, these quantities are fed in to the sum- 
over-states formulas of Orr and Ward(2& to obtain the 

(3) 

correlated values of Xthg- More details about the pro- 
cedural aspects of various CI approaches used by us can 
be found in our earlier works £L2iL22i21 



The total magnitude of component of the longitudinal 
component of THG susceptibility (\Xxxxx(—3v;v,v,l))\), 
as a function of the incident photon energy E, for fifty 
unit cell oligomers of PDPA and irans-polyacetylene is 
plotted in Fig. |2j From the figure it is clear that the main 
peak of the spectrum for both these materials is the first 
one corresponding to the lowest optically allowed state 
1B U , and is located at Huj = E g /3, where E g is the value 
of the optical gap defined as E g = E(1B U ) — E(lA g ). 
From the figure it is obvious that: (a) qualitatively speak- 
ing, longitudinal THG spectra for the both the polymers 
are similar (b) from a quantitative point of view, as ex- 
pected, the resonant features of the PDPA spectrum are 
red shifted as compared to the irans-polyacetylene ones 
(c) in the frequency region of interest, both the reso- 
nant, and the off-resonant, nonlinear optical responses 
of PDPA are much more intense as compared to trans- 
poly acetylene. Since the 1B U state of PDPA's in Hiickel 
model originates from the excitation HOMO — > LU MO, 
both of which are chain-based orbitals, one can safely 
conclude that as far as the longitudinal nonlinear-optical 
response computed at the Hiickel level is concerned, 
oligo-PDPA's behave as if they were smaller band-gap 
polyenes. 



2. Transverse Component 



III. RESULTS AND DISCUSSION 

(3) 

In this section we present our results on Xthg usm g 
various approaches. First we present and briefly discuss 
the results computed using the Hiickel model, followed 
by the results computed using various CI approaches. 



A. Independent-electron theory 

1. Longitudinal Component 

Here we present and discuss the longitudinal THG 
spectrum for the case of oligomer PDPA-50, which has 
700 sites leading to 350 occupied and the same number 
of virtual orbitals. Even within the one-electron theory, 
calculation of x^ ' s with such a large number of orbitals 
is intractable. Therefore, to perform such calculations 
on larger oligomers, it is mandatory that some of the 
orbitals be discarded. To this end we adopted the crite- 
rion to discard those occupied and virtual orbitals which 
are more than 3.0 eV away from the Fermi energy. The 
choice of 3.0 eV is based upon the magnitude of the op- 
tical gap and the location of various features in the x' 3 ^ 
spectra of smaller oligomers. Since, on PDPA-30 it was 
also possible to perform exact calculations, we were able 
to compare our approximate results with the exact ones 
for that oligomer. We found that the approximate results 
on longer oligomers are quite accurate. 



The transverse THG spectrum of oligo-PDPA's satu- 
rates very rapidly with size because with increasing con- 
jugation length, the irans-stilbene-like structure of the 
oligomer in the transverse direction is unchanged. With 
increasing conjugation length, only the overall intensity 
in the transverse direction is expected to increase, with 
little change in the qualitative feature. Therefore, here 
we present the exact results on \Xyyyy(—3w;u),U),U))\ for 
the oligomer PDPA-30 in Fig. 03 It is clear from the 
figure that there is only one resonant feature in the spec- 
trum and it is located close to 1.7 eV and its intensity 
is about one-fifth that of the largest resonant intensity 
in the longitudinal THG spectrum (cf. Fig. |2j. Even 
for smaller oligomers such as PDPA-10 and PDPA-20, 
the location of the resonance is the same pointing to 
its possible origins in the phenyl-based orbitals. There- 
fore, to simplify our task of locating these orbitals, we 
investigate the THG spectrum of PDPA-10 for which the 
main peak occurs at 1.71 eV, corresponding to an A g 
type state, obtained from the ground state by single- 
excitations H — ► L + 39 and H — 39 — > L. To in- 
vestigate the nature of the orbital corresponding to the 
L + 39-th one-electron state, we calculated the contri- 
bution of the charge density centered on the backbone 
carbon atoms, to its total normalization. The contribu- 
tion was computed to be 0.05 which indicates that the 
orbital in question is predominantly centered on the side 
phenyl rings. Further investigation of the orbital coef- 
ficient reveals that the orbitals in question, are derived 
from the phenyl-based delocalized (d*) virtual orbitals, 
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with no contribution from the localized orbitals (7*-type) 
of the phenyl rings. Similarly, owing to the particle-hole 
symmetry, H — 39-th orbital is derived from the d-type 
occupied orbitals of the phenyl rings. The investigation 
of longer oligomers yields identical results. Thus it is 
clear that the transverse THG susceptibility of the sys- 
tem owes its origins to phenyl-based delocalized (d/d*) 
levels. This is an important point which also helps us 
perform effective correlated calculations of this compo- 
nent, presented in the next section. 

B. Correlated-electron theory 

Since the smaller energy gaps obtained with the 
screened parameters in our earlier works were found 
to be in much better agreement with the experiments, 
we will present our main results based upon screened- 
parameter-based calculations. However, when we com- 
pare the PDPA nonlinear optical spectra with those of 
polyenes, we will use the standard parameters because 
screened parameters are not valid for polyenes. 

1. Longitudinal Component 

Now we present correlated calculations for the lon- 
gitudinal component of the third-harmonic generation 
susceptibility Xxxxx{— 3w; u>, u>, w) for oligo-PDPA's, per- 
formed using the QCI method. Given the large number 
of electrons in these systems, QCI method is not feasible 
for them if all the orbitals of the system are retained in 
the calculations. Since the longitudinal nonlinear opti- 
cal properties are determined by low-lying excited states 
of the system, in the limited CI calculations we decided 
to include the orbitals closest to the Fermi level. There- 
fore, for PDPA-n we included n occupied, and n vir- 
tual orbitals closest to the Fermi level in the QCI cal- 
culations. Remaining occupied orbitals were frozen and 
virtual orbitals were deleted as explained in section [H] 
Thus, the computational effort associated with the QCI 
calculations on PDPA-n is same as that needed for a 
polyene with n double bonds. Although for PDPA-10, it 
leads to Hilbert space dimensions in excess of one million, 
however, using the methodology reported in our earlier 
work o 1 *?' 2 ^ , we were able to obtain low-lying excited states 
of such systems. 

The longitudinal THG spectra (\xxxxx(— 3ui;ui,ui,ui)\) 
for PDPA-5 and PDPA-10 employing the screened pa- 
rameters and the QCI method are presented in Figs^Ja) 
and 0Jb) , respectively. From the figures it is obvious 
that the spectra for PDPA-5 and PDPA-10 containing 
four main features each are qualitatively quite similar, 
suggesting the possibility that an oligomer as small as 
PDPA-5 may possess essential features of bulk PDPA. 
It is clear from Fig. 21 that: (a) The intensity of the 
nonlinear response increases with the increasing length 
of the oligomer, and (b) the peaks of the oligomers are 



redshifted with the increasing conjugation length. Next 
we examine these peaks in detail. 

The properties of excited states contributing to the 
longitudinal THG spectra of PDPA-5 and PDPA-10 are 
presented in table Q] From the table, and Fig. 31 it is 
clear that there are three £?„-type states, 1B U , nB u , and 
kB u along with two A 9 -type states, 2A g and mA g which 
contribute, respectively, to the three-photon and two- 
photon resonances in the spectra. Upon comparing the 
wave functions of various excited states of PDPA-5, with 
those of PDPA-10, we infer that qualitatively the wave 
functions are very similar. For each oligomer, 1B U state, 
which constitutes the first peak of the spectrum consists 
mainly of across-the-gap, single excitation. The next 
peak corresponds to the nB u state which, for both the 
oligomers, also consists mainly of higher energy single ex- 
citations. The third three-photon resonance correspond- 
ing to a still higher energy state labeled kB u , is a mix- 
ture of both singly- and doubly-excited configurations. 
For PDPA-10, this state is dominated by a high-energy 
single excitation with important contributions from dou- 
ble excitations, while in PDPA-5 the situation is reverse 
with main contributions coming from the double excita- 
tions. Next, when we examine the wave functions of the 
states 2A g and mA g contributing to the two-photon fea- 
tures, we find that configurations contributing to them 
are essentially identical. However, because of the orthog- 
onality constraint, the relative signs of the coefficients of 
these configurations are opposite for the two states. One 
strange aspect of the THG spectrum presented in Fig. 
His the coincidence of the two-photon and three-photon 
resonances corresponding to states nB u and 2A g , respec- 
tively, for both PDPA-5 and PDPA-10. This is because 
for both the oligomers, E(nB u )/3 « E(2A g )/2. One 
wonders whether this is an artifact of using a truncated 
orbital set in the QCI calculations, or a genuine effect. In 
order to perform correlated calculations of the transverse 
THG spectrum of PDPA-5, we used an extended orbital 
set and performed large-scale MRSDCI calculations pre- 
sented later in Fig. From there it is clear that although 
both 2A g and nB u make important contributions to the 
spectrum, however, their peaks are well separated, and 
the intensity of nB u is much higher than that of 2A g , 
leading us to conclude that the coincidence of 2A g /nB u 
peaks in the QCI/THG spectrum is an artifact. 

Next we compare the longitudinal THG spectrum of 
PDPA-10 with that of a polyene consisting of ten double 
bonds (PA-10), computed using the QCI method, and 
the standard parameters in the P-P-P Hamiltonian. The 
spectra are presented in Fig. while the properties of 
the excited states contributing to the spectra are illus- 
trated in table |Ul It is clear from that the intensity 
of the PDPA-10 spectrum is generally larger and all its 
major peaks are significantly redshifted as compared to 
the irans-polyacetylene oligomer, in agreement with the 
Hiickel model results. Moreover, the THG spectrum of 
PDPA-10 computed with the standard parameters, still 
has four main peak corresponding to the states 1B U , 
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nB u , kB u , and mA g , in agreement with the screened- 
parameter results. The main difference in the THG spec- 
trum of PDPA-10 computed with standard parameters, 
as compared to the one computed with the screened pa- 
rameters, is that there is no contribution to it from the 
2A g state. When we compare this aspect of PDPA-10 
standard parameter THG spectrum with that of PA- 
10, we find that the PA-10 spectrum also does not ex- 
hibit any feature corresponding to 2A g . Moreover, PA- 
10 THG spectrum has only three peaks corresponding 
to 1B U , nB u , and mA g , with no peak corresponding to 
a fc-B„-type state. If we compare the wave functions of 
these three states of PA-10 with those of PDPA-10 (tar 
blelnj, we find that the wave functions of 1B U and mA g 
states are very similar to each other for the two materials. 
As far as the nB u state is concerned, its wave function 
for PDPA-10 is composed predominantly of single exci- 
tations, while for PA-10 it exhibits large configuration 
mixing, with significant contribution also coming from 
the double excitations. Moreover, in PDPA-10, the nB u 
state is very close energetically to the mA g state (slightly 
below it), while in PA-10, it is much higher in energy 
than the nB u state. Note that the nB u state of PDPA- 
10 computed with the screened parameters is slightly 
above the mA g state (table Pi . Therefore, the location 
of the nB u state is an interesting difference between PA- 
10 and PDPA-10 which these standard-parameter-based 
calculations suggest, and it merits further investigation 
in possibly more extensive calculations. Thus, experi- 
mental investigations of the longitudinal THG spectrum 
of PDPA-10 could possibly shed some light on the follow- 
ing aspects: (a) whether or not the 2A g state makes an 
important contribution to the spectrum indicating as to 
what range of Coulomb parameters (screened/standard) 
are valid for the material, (b) location of nB u vis-a-vis 
mA g , and (c) whether or not there is a peak correspond- 
ing to the kB u state in PDPA's. 



2. Transverse Component 
Performing accurate correlated calculations of 

(3) 

Xyyyy{—3u>;uj,uj,uj) for oligo-PDPA's is an extremely 
difficult task. The reason being that the many-body 
A g -type states which contribute to the peaks in this 
component of the susceptibility are very high in energy, 
and hence are more difficult to compute by many-body 
methods, as compared to the low-lying excited states 
contributing to the longitudinal spectra. The very high 
excitation energies of these states are due to the fact 
that these A g -type states are predominantly composed 
of excited configurations involving high-energy delocal- 
ized orbitals originating from the phenyl substituents. 
However, during our independent-electron study, we 
concluded that transverse susceptibilities exhibit rapid 
saturation with the conjugation length. Therefore, 
we restrict our study of the transverse components of 
the nonlinear susceptibilities to PDPA-5. We perform 



calculations with the screened parameters and the 
MRSDCI method employing thirty orbitals in all, of 
which fifteen were occupied orbitals, and the remaining 
fifteen were the virtual ones. Rest of the occupied 
orbitals were frozen. Of the thirty orbitals, ten were 
the orbitals closest to the Fermi level which were also 
used in the QCI calculations. Remaining twenty orbitals 
were the d/d* -type phenylene-ring-based orbitals closest 
to the Fermi level. In the MRSDCI calculations, no 
l/l*-type phenylene orbitals were used because they do 
not make any significant contribution to the transverse 
nonlinear spectra in the Hiickel model calculations. 
In the MRSDCI calculations, we used 25 reference 
configurations for the A g -type states, and 24 for the 
Su-type states leading to CI matrices of dimensions 
close to half-a-million both for the A g and B u manifolds. 

The transverse THG spectrum of PDPA-5, along with 
the THG spectrum of irtms-stilbene (labeled PPV-2) 
oriented along the y-axis computed using the QCI ap- 
proach and screened parameters, are presented in Fig. 
El The properties of the excited states contributing to 
these spectra are displayed in table IIIII Before dis- 
cussing the many-particle wave functions of various ex- 
cited states, it is important to specify that for PDPA-5, 
orbitals H — 17 /L + 17 and H — 19 /L + 19 appearing in 
various states are the substituent phenyl-ring-based d/d* 
orbitals making important contributions to the trans- 
verse THG spectra computed using the Hiickel model 
(cf. sec. IIII A2|) . It is obvious from Fig. EJthat: (a) 
For PDPA-5 four states labeled n y B u , k y B u , j y A g , and 
rriyA g (subscript y implies y, i.e., transverse direction) 
contribute to the spectrum, while (b) for PPV-2 three 
states 1B U , jA g , and mA g make main contributions. 
Furthermore, the main contributions to the intensities 
of the THG spectra of both the substances are from two 
Ag-type peaks: j y A g and m y A g for PDPA-5, and jA g 
and mA g for PPV-2. For PDPA, we recall that at the 
Hiickel model level, the transverse THG spectrum had 
only one peak due to closely-spaced A g -type states cor- 
responding to a single-particle excitation from the or- 
bitals near the Fermi level (H/L) to the phenyl-based 
d/d* orbitals. Upon examining the many-particle wave 
functions of various excited states of PDPA-5, we con- 
clude that even at the many-body level the excited states 
contributing to the transverse THG spectra originate 
from the same set of single excitations involving H/L 
and phenyl-based d/d* orbitals except for the k y B u s,tsXe. 
However, the charge-density analysis of the orbitals in- 
volved in the k y B u state reveals that the orbitals in ques- 
tion (H — 3/L + 3, H — 4/L + 4) have significant charge 
density on the substituent phenyls. Thus, we can safely 
conclude that in the transverse THG spectrum of PDPA- 
5, the substituent-phenyl-based orbitals play a very im- 
portant role. On comparing the many-particle wave func- 
tions of the Ag-type states contributing to the THG spec- 
tra of the two systems, we find that unlike the case of 
PDPA-5, in PPV-2 these states (kA g and mA g ) have im- 
portant contributions from the two-particle excitations 
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as well. However, the contributions of the single-particle 
excitations to these states in both PDPA-5 and PPV-2 
have similar character — they involve H/L orbitals and 
high-lying orbitals localized predominantly in the phenyl 
rings. Upon further investigation of these high-lying or- 
bitals involved in single excitations we conclude that they 
are of d/d*-type phenyl-based orbitals. Therefore, we 
conclude that the mechanism of the transverse THG in 
PDPA's has noticeable similarities to that of THG in 
irans-stilbene. 

Finally, in Fig. we compare the transverse THG 
spectrum of PDPA-5 with its longitudinal one, both com- 
puted using the MRSDCI procedure, and the screened 
parameters. It is clear from the figure that the most in- 
tense peak in the transverse spectrum is about half in in- 
tensity as compared to the strongest peak in the longitu- 
dinal one, and is located in the higher energy range where 
there is no competing longitudinal response. Thus, our 
PDPA-5 calculations suggest that the transverse THG re- 
sponse of these materials, compared to their longitudinal 
response, is significant, as well as clearly distinguishable 
from it, in an experiment. However, it remains to be 
verified whether the results obtained for small oligomers 
will also be valid for the bulk PDPA's. 



IV. CONCLUSIONS 

Our aim behind undertaking the present theoretical 
study of the nonlinear optical properties of the novel 
polymer PDPA was, not only to calculate their THG 
spectra, but also to understand the underlying mecha- 
nism in a way similar to what has been possible for sim- 
pler polymers such as irens-polyacetylenep^ PPP, and 
ppy | 9 i 24 This understanding has chiefly arisen from the 
essential state mechanism alluded to earlier in section [IJ 
which aims to explain the nonlinear optical properties of 
chain-like ^ and more recently phenylene-based conju- 
gated polymers^ in terms of a small number of excited 
states. For PDPA's, whose structure has ingredients in 



common with both the chain-like as well as the phenyl- 
based polymers, our calculations suggest that their non- 
linear optical properties can also be understood along the 
same lines as the essential state mechanism. However, 
because of the anisotropic nonlinear optical response of 
PDPA's stemming from their structural anisotropy, the 
essential states contributing to various components of 
susceptibility are different. Our independent-electron 
as well as correlated calculations suggest that the THG 
spectrum of PDPA's is distributed over two natural en- 
ergy scales: (a) low-energy region of the spectrum has 
significant intensity mainly for the longitudinal compo- 
nent, while (b) high-energy region of the spectrum has 
intensity mostly for the transverse component. The lon- 
gitudinal THG spectrum bears striking qualitative re- 
semblance to the THG spectrum of polyenes with simi- 
lar essential states contributing to the intensities therein. 
The transverse THG spectrum on the other hand, has re- 
semblance to the THG spectrum of irans-stilbene, with 
again, similar type of essential states contributing to their 
intensities. The essential states which contribute to the 
longitudinal spectra of these polymers are quite distinct 
from the states contributing to their transverse spectra. 
The states contributing to the longitudinal spectra can 
be described in terms of chain-based orbitals close to the 
Fermi level, however, to describe the transverse spectra, 
one needs to take into account the excitations involving 
phenyl-based orbitals away from the Fermi level. There- 
fore, it will be of considerable interest if the theoretical 
predictions presented here can be tested in experiments. 
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Figure 1: The unit cell of PDPA. The phenyl rings are rotated 
with respect to the y-axis, which is transverse to the axis of 
the polyene backbone (a>axis) 
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Figure 2: Comparison of \xxxxx(—3cu; uj, u>, lo)\ ofPDPA (solid 
line) and trans-polyacetylene (dotted line) oligomers con- 
taining fifty unit cells, computed using the Hiickel model. 
Linewidth of 0.05 eV was assumed for all energy levels. 
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Figure 3: \Xyyyy(~ 3<^; uj, w)\ of PDPA-30 plotted as a func- 
tion of incident photon energy. Linewidth of 0.05 eV was as- 
sumed for all energy levels. 




Figure 4: The magnitudes of the longitudinal component of 
the third-harmonic generation spectra (|XsmLs( — 3lu; lo, uj, uj)\) 
of oligo-PDPA's computed using the QCI approach and 
the screened parameters for: (a) PDPA-5 (b) PDPA-10. A 
linewidth of 0.05 eV was assumed for all the levels. 



Table I: Excited states contributing to the longitudinal THG 
spectrum of PDPA-5 and PDPA-10 as obtained in QCI cal- 
culations. Under the heading wave function, we list the most 
important configurations contributing to the many-body wave 
function of the state concerned, along with their coefficients, 
consistent with our convention--. 





PDPA-5 


PDPA-10 


State 


Energy (eV) 


Wave Function 


Energy (eV) 


Wave Function 


2A g 


2.86 


|ff -> L + 1) +c.c. (0.57) 
|ff -» L;H -> L)(0.48) 


2.08 


|ff -> L + 1) + cc. (0.49) 
|ff -> L;ff -> L)(0.50) 


mA g 


3.98 


|ff -fL + l)+c.c.(0.37) 
|ff -> L;H -> L)(0.74) 


2.82 (3A g ) 
2.88 (4A 9 ) 


|ff-> L + l) + c.c.(0.68) 
|ff->L;ff->L)(0.37) 


1B U 


2.59 


|ff -> L)(0.99) 


1.98 


|ff -> L)(0.96) 


nB u 


4.21 


|ff- 1 -> L+ 1)(0.94) 

|ff -> L; H -> L + 1) + c.c.(0.14) 


3.08 


|ff -1-+L + 1) (0.72) 
|ff -> L + 2) +c.c.(0.38) 


kB u 


5.58 


|ff — > L; ff — > L + 1} + cc. (0.57) 
IH-l-fX + ljH-l-f^H- c.c.(0.25) 
|ff-l -»L+1)(0.19) 


3.96 


|ff-2^i + 2)(0.59) 

|ff -> L; H -> L + 1) + c.c.(0.31) 

\H - 1 -> L + 1;H - 1 -> L) + c.c.(Q.23) 
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Figure 5: The magnitude of the third-harmonic generation 
spectrum (\xxxxx{~ 3lu; lo, lj, u))\) of: (a) PDPA-10 (top) and 
(b) ten unit oligomer of trans-polyacetylene (bottom) com- 
puted using the QCI approach and the standard parameters. 
A linewidth of 0.05 eV was assumed for all the levels. 



Table II: Comparison of excited states of ten unit polyene 
(PA-10) and PDPA-10 computed using the QCI method and 
the standard parameters in the P-P-P Hamiltonian. Rest of 
the information is the same as given in the caption of table |l| 





PDPA-10 


PA-10 


State 


Energy (eV) 


Wave Function 


Energy (eV) 


Wave Function 


IB U 


2.91 


\H -> L>(0.90) 
\H-1 -> L+ 1)(0.32) 


3.53 


|ff -> L)(0.82) 
|ff-l -> L + l) (0.28) 


mA g 


5.05 


\H ->L;H -> L)(0.47) 
|ff--»L;ff-l-» L + l) (0.45) 
\H -> L + l) +c.c.(0.33) 


5.49 


|ff->L + l) + c.c.(0.42) 

|ff -» L;H -> L)(0.36) 

|ff -> L;ff - 1 -> L+ 1)(0.31) 


nB u 


4.90 


1 -> L + 1)(0.62) 
|ff ->-L + 2)+ c.c. (0.29) 


7.31 


|ff -> L; H -»• f, + 1) + c-c. (0.28) 
|ff -> L; ff - 1 -» L + 2) + c.c.(0.27) 
|ff-l-> L + l) (0.23) 


kB u 


6.65 


|ff -»• L;H -> L + 1) + e.c.(0.36) 
|ff- 1 -f L + l) (0.29) 
|ff -2^L + 2>(0.20) 
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Figure 6: Comparison of transverse third-harmonic genera- 
tion spectra of PDPA-5 with that of trans-stilbene (PPV-2) 
oriented along the y-axis. Both the spectra were computed us- 
ing the screened parameters. For PDPA-5, MRSDCI method 
was used, while for PPV-2, QCI method was used. A linewidth 
of 0.05 eV was assumed for all the levels. 



Table III: Nature of excited states contributing to the trans- 
verse THG spectrum of PDPA-5 computed by MRSDCI 
method. Rest of the information is same as given in the cap- 
tion of table |l| 



PDPA-5 


PPV-2 


State 


Energy (eV) 


Wave Function 


State 


Energy (eV) 


Wave Function 


jyAg 


4.39 


\H -> L + 17) + c.c (0.43) 
\H -> L + 19) + c.c (0.43) 


jA g 


4.45 


\H -> L + 3) + c.c (0.43) 

L-H ^ L} + c.c (0.42) 


TUyAg 


4.52 


\H -> L + 17) + c.c (0.27) 
\H -> L + 19) + c.c (0.38) 
\H - 2 -> L + 1) + c.c (0.30) 


mA g 


5.72 


\H - 1 -> L + 1;H -> L)(0.29) 
|Jf-2^L + 2;//^ L)(0.29) 
|tf ->Z + 3) + c.c. (0.26) 
|.ff -> L + 5) +c.c.(0.25) 


Fly 


4.96 


\H -> L + 19) +c.c.(0.60) 


IB U 


4.35 


|H -> L)(0.90) 


kyBu 


5.61 


|#-3^L + 3)(0.56) 
|#-4^L + 4)(0.49) 
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Figure 7: Comparison of longitudinal (solid lines ) and the 
transverse (dotted lines) third-harmonic generation spectra 
of PDPA-5 computed using the screened parameters and the 
MRSDCI method. A linewidth of 0.05 eV was assumed for all 
the levels. 



